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ABSTRACT 

There is increasing interest in extended source effects in microlensing events, as probes 
of the unresolved sources. Previous work has either presumed a uniform source, or else 
required an approximate or numerical treatment of the amplification function averaged 
over the source disk. In this paper, I present analytic expressions for the angle-averaged 
amplification functions for the rotationally-symmetric intensity and polarization cases. 

These integrals will allow us to use the technology of inverse problems to study 
the source limb-darkening and limb-polarization functions. 

Key vifords: gravitational lensing - methods: data analysis 



1 INTRODUCTION 

The most common interest in microlensing events is as a 
probe of the lensing object itself. Recently, however, there 



otherwise unresolved sources ( Valls-Gabaud 1998 




Mao and 


Witt 1998 




Gaudi and Gould 1999]). When the lens transits 



the source (or nearly so), it breaks any rotational symme- 
try, and this giv es us access to the sur f ace brightness, as wel l 



as polarization (Simmons ct al. 1995 



immons et al. 1991: 



and chromaticity ( Valls-Gabaud 1995{ Valls-Gabaud 199J) 
information. Previous work on extended source effects has 
concentrated on the forward problem and generally either 
performed the required calculations numerically rather than 
analytically, or used an approximate form of the amplifica- 
tion function. Also, much of the work on source effects has 
relied on the high amplification provided by binary lens cau- 
sic crossings, rather than the amplification of a single lens. 

The gravitational lensing forward problem - that of 
predicting centroid motion and magnification for a given 
set of source parameters - is relatively easy. The problem 
is also, however, typically poorly conditioned, in the sense 
that there will be a broad range of limb-darkening or limb- 
polarization functions which could plausibly correspond to 
the observed signal in a microlensing event. This means that 
a parameter-fitting approach to recovering these functions is 
very dangerous. 

We can make progress by expressing the problem explic- 
itly as an inverse problem (IP), and using the technology of 
IP methods to analyse precisely what information can be 
recovered for a given set of observations. 

Th at is the subject of, and m otivation for, a forthcoming 
paper ( Sray and Coleman 200C ); here I am concerned with 
identifying the angle- averaged amplification functions as IP 



kernels, and obtaining analytic expressions for them. As well 
as facilitating the IP analysis of the problem, these analytic 
kernels can help in the treatment of the forward problem, 
since they can be evaluated more efficiently than by a nu- 
merical integration, and with high accuracy over their entire 
domain. 

In Sect. ^, I define the amplification functions as IP 
kernels, in Sect. ^, I integrate them about the source's centre, 
and in Sect. 0, I present the results of that integration. 



2 AMPLIFICATION FUNCTIONS AS 
INVERSE PROBLEM KERNELS 

The geometry of a microlensing event is as shown in Fig. [l|. 



The gravit ational lens amplification function is (Schnei 



der et al. 1992) 



(1) 



where 



1/2 



C^^ = r'^ + s'^ — 2rscos(x — 0)- 



Denote the intensity by I{r) and the Stokes parameter 
by Q{r,x) = — P(r)cos2x, where P{r) is the polarization 
of the stellar surface, and we are assuming that the sur- 
face brightness is rotationally symmetric. In the case of a 
microlensing event, we cannot resolve details of the lensed 
source, and must therefore measure integrals over the source 
surface. We immediately obtain 



I{s{t),^{t)) = 



I{r)Ai{r- 



(2) 
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Figure 1. Geometry of a lensing event. The projected path of 
the lens has impact parameter /, and the path is parameterised 
by polar coordinates s(t) and (pit), relative to the centre of the 
source, projected into the lens plane. Any point in that plane can 
be given in polar coordinates r and Xi and this point is a distance f 
from the centre of the lens. All dimensions are normalised to the 
Einstein radius in the source plane. The angles <j> and x are taken 
with respect to the line joining the source to the lens' point of 
closest approach. 



Qis{t),<t>{t)) 



P{r)AQ{r; s, <fi) dr 



where the amplification kernels are 



Ai = r A{r,x;s,(t))dx 



Aq = -r / cos 2xA(r, x; s, </>) dx- 
Jo 



(3) 

(4) 
(5) 



Note that the kernel Ai is a factor 2Trr times the angular av- 
erage of the amplification function, and the functions 7(s, 0) 
and Q{s,<f)) have the dimensions of fiux rather than inten- 
sity- 

Equations (H) and (13) are in the form of an inverse 
problem. We address the inverse p roblem in a forthcom- 
ing paper (Gray and Coleman 2000). The evaluation of the 
integrals Ai^q is rather hard, and I describe it in this paper. 



3 INTEGRATION OF THE AMPLIFICATION 
FUNCTIONS 

Write 

z — exp i{x — 4') , dz — izdx, (6) 
so that 

a2 2,2 / I 1\ 

Q = r + s — rs I z H I . 

Define 



ai = 
and 



r^ + s^ 
2rs ' 



ai 



4 + -I- 



2rs 

r / s ,r > s 
s/r ,r<s 



X2 — £12 + \/a| 
Defining Xi = Xjxi, we have 



(7) 

(8) 
(9) 



Figure 2. The contour for the integral in Eqn. (p^, showing the 
cuts between the poles at a;2, xi, x\ and X2- 



Xi -\- Xi — 2£lj , X-iXi — 1. 

It is easy to show that 

< a;2 < a;i < 1 < a::i < a;2. 

Rewriting the expression for ^, we obtain 



4 



1/2 



{z — X2)(z — X2) 



_ x \ 1/2 



{z - Xl)(z — Xl) 

We now evaluate the integral 



(10) 



(11) 

The integrand has 



where C is the contour shown in Fig. 
poles at z — xi,2, z — xi,2 and z — 0. 

The contour encloses a single singularity at 2: = 0. We 
have ^(2: = 0) = 1, so that 

/ = 27ri X Res(7(2 = 0)) = 47ri. (12) 

For contour 1, substitute z — e'^ , for ip G [0,27r]. Then 
/f-*^' = i J^^ {$, + 1/^) dip, and the substitution tp = x ~ 
produces 

/W=2i/ A{ax,(P))dx=-Ai{r;s,(P). (13) 

J <f) 

Contours 5 and 6 cancel, and with the substitution z — 
X2 + pe'*, it is clear that j'^-* ^ as p — > 0. 

Now turn to contours 2 and 4. By substituting z = X2 + 
a, a = ^ {xi —X2), into substituting z = X2 + cre'^^ , 
a — (xi -~X2), into /'^', and noting that < a < a;i — 5^2 
(so that \a\ < \x2 ^ X2\, \x2 — xi\, and \x2 ^xi\, so that the 
phases of the corresponding square-rooted factors in f'-^-' are 
unaffected by the factor of e''^'^), we can see that 

7(^' (14) 

Now substituting z — x {x real) directly into Eqn. (|l^), we 
obtain 



where 



h = 
h = 



{X2 — x){x — X2) 

X \ {xi — x){xi — x) 

1 / {X2 — x){x — X2) 



1/2 



da; 



-1/2 



da;. 



(15) 

(16) 
(17) 



(a;i — a;)(a;i — 3 

© 0000 RAS, MNRAS 000, 000-000 



Angle-averaged microlensing amplification functions 3 



Using the notation of ( [Carlson 1988 ) 

p n 

[pi,...,p„] = j Wia^ + htf^^Ut, 



we may rewrite these as 

h = [+1,-1,-1, +1,-2] 
h = [-l,+l,+l,-l,-2] 



(18) 
(19) 



with fli — {—X2,xi,xi,X2,Q) and bi — (+1, —1, —1, —1, +1) 
in both cases. We will evaluate these integrals when we re- 
turn to them below, in Sect. W. 

Thus, collecting Eqns. (^, (Q, (Q and ([s]), we ob- 
tain the angle- averaged amplification function as 



Ai{r;s) =r(27r + Ji ~ h 



(20) 



This has been confirmed by direct numerical integration of 
the integrand. 

Turning now to Eqn. (|^, we may again substitute z = 
&q)i(X ~ 0)i and obtain 

r, 1 / 2 i2<ft , -2 -i2rf>\ 

cos lx= -^\z e + z e j . 
Now evaluate the integral 







i2<h , -3 - 

ze + z e 



dz, 



(21) 



with the same contour as above, and with ^ as in Eqn. (p^. 
This has a third-order pole at 2: = 0, which means that the 
residue is 



1 d^ 
Res(Q) = a„i ^ 



with ai as defined in Eqn. 
Q = 27rie~'^'*(ai — 02)^^. 



= (ai-a2)" (22) 
above. Thus 



(23) 



Much of the calculation goes through as before. Substi- 
tuting z = e^^ , we obtain, 



Ai - 

= AQ{r;s,(i>). 



(24) 



Contours 5 and 6 cancel, and contour 3 makes zero contri- 
bution in the p — * limit. Similarly, 

q(2) ^ q{4) ^ ie''"^(-Qi + Q2) + ie-'^'^(-g3 + Q4), (25) 
where 



Qi = [+i,-i,-i,+i,+2l 

Q2 = [-l,+l,+l,-l,+2] 

Q3 = [+1,-1,-1, +1,-6] 

Q4 = [-l,+l,+l,-l,-6] 



(26) 
(27) 
(28) 
(29) 



with Oi and bi as above. 

Assembling equations (E3|), (e4) and (Ea), we obtain 



AQ{r-s,(t>) 



[e'^^-Oi + Q2) + e-'''^(-Q3 + Q4) 
— 7re~'^'*(ai — 02)^] • 



However, this integral should be real, so the imaginary part 
must be zero: 

— Qi + Q2 + Q3 — Q4 + 7i"(ai — 02)^ = 0. 




0.8 1.0 



Figure 3. The amplification kernel A/(r; s), plotted as a function 
of r for a selection of values of the distance s, and for an impact 
parameter I = 0.1 (see Fig. There is a singularity along the 
line r = s, where the integration includes the point C = 0- 



Thus, the final expression for the angle-averaged polariza- 
tion amplification function is 



AQ{r;s, 



= — r cos 2(f>{Qi — ' 



(30) 



1 (Q1-Q2) 



This also has been confirmed by numerical integration. 

Eqn. ( pc| ) and Eqn. (^) are the principal results of this 
paper. 

In Fig. ^ I show the amplification function Ai{r;s), 
with a singularity along the line r = s. Note particularly 
the breadth of the kernel after the singularity: although the 
angle- averaged amplification is close to 1 away from a well- 
defined peak at r — s, the factor of r in Eqn. (^) means 
that there are contributions to I{t) in Eqn. (^) from a broad 
range of 7(r), a situation especially severe for cases where the 
source function I{r) extends significantly beyond r — 1 {le, 
those cases where J?« > Re)- The breadth and asymmetry 
of the kernel is what makes the recovery of the source func- 
tion so problematic. Similarly, Fig. ^ shows the amplification 
kernel AQ{r; s, (j)). This has the same factor of r as Ai{r; s), 
but because the underlying angle- averaged function is close 
to zero away from r = s, this has a less damaging effect, 
so that the polarization signal should be easier to recover 
(which is fortunate, since that signal is so much harder to 
detect than the intensity signal) . The kernel is still, however, 
both broad and asymmetric. 



4 AMPLIFICATION FUNCTIONS AS 
ELLIPTIC INTEGRALS 

The integrals defined by equations (^) and (^|) , and equa- 
tions (Eg) to (E9), are elliptic integrals. Carlson (1988) pro- 



vides a set of recurrence relations to reduce such integrals 
to a small set of elementary integrals, which are in turn ex- 
pressed in terms of a set of functions Rj, Rf and Rd which 
are more symmetrical alternatives to the traditional Legen- 
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Figure 4. The amplification kernel AQ(r;s,(l>), plotted as in 
Fig. ^; the parameter <p(s) is that corresponding to a straight- 
line transit with impact parameter i = 0.1. 



dre elliptic integrals, and which can be evaluated efficiently, 
and to hig h accuracy, across their entire domain. A suitable 
algorithm ( [Gray 200(]| ) can do this (algebraically exhausting) 
reduction mechanically, to obtain 



Ai{r-s) 



2-K 



h 



+ 3-^2(1 + \\)^X\Xl 



Rj2 

X2 



(31) 



2 ^' I^ A'(k) 



^/xix^ 



+ 4-!-^^V^(n2-ni), (32) 

X2 



— r cos 2(j) 



(^1 - (^2 

_ qpRp + qpRp + qjR.J2 
2(a;ia;2)3/2 

_ qeEjn) + qkK{n) + 9^112 
2(ccia;2)3/2 

Here 

A, . ^ii^>0, 

a;i — X2 

A2 = ^1^>0, 
2:2 — a;i 

K = A1A2 (0 < < 1); 
the coefficients qa are 

(xi + x2){xxX2 + l)(a;2 — xxY 



(33) 
(34) 



go = 



qp = 



3(0:1X2 — 1) 

{xiX2 — 1)(2;2 — XiY 
X2 

„2 „ n3 



(xj-l) (3:2- 3:1)' 
3X2 



(7e = 

qk = 
q-TT = 



3<3i5 / , \ / 2 2 1 X 

5- = + X2)[XiX2 - 1), 



3gd , , 3g.j 
^ + '^^ + — 

(—2:2 + 3a:ia:2 + xi + X2){xiX2 — 1), 

3g.j _ (3:2 - l)(a:i3:2 - l)(a:2 - a:i)^ ^ 
Ai 3:2 



we have used the abbreviations 
Rp = i?F(0, 1 - K^ 1), 

i?D = i?D(0, l-K^l), 

i?,7> = i^J(0,l-«:^l,l + AO, 



n 



and n(0 = vr/2, n, k), K{k) and E{k) are the Leg endre forms 
of the complete elliptic integ rals, as defined in ( Sradshte^ 



and Ryzhik 1994, 8.111), or (Abramowitz and Stcgun 1965 



17.2) (with k = sin a). These are related to the Carlson 
forms of the integrals by 

7?j(0,1-k',1;1 + A) = hK(K)-n{7v/2,-X,K)] 

A 

Rd{0,1~k,1) = -^[K{ti)-E{n)] 
i?F(0, l-«:^l) = K{k). 



4.1 Singularies and asymptotic behaviour of Aj 

and Aq 

We c an write Eqn. (p^) in terms of H euman's Lambda func- 
tion (Abramowitz and Stcgun 1965, 17.4.39) as follows: 



2 ^''^ K{k) 



^X\X2 

+ 27r 



Ao(arcsin(l + Ai)"^/^ \ a) 



(Ai ^ A2) 



(35) 



The advantage of this is that we can now easily isolate the 
singularity at r = s, where we have a:i = 3^1 = 1, 3:2 > 1, and 
thus K = \. The Aq function has no singularities, and the 
coefficient of K{k) i s finite there, so the only sing ularity is 
at K{k = 1) where (Abramowitz and Stcgun 1965, 17.3.26) 



Thus the leading order term in Eqn. (B2|) at r = s is 



Ai(r ~ s; s) 



2-K 



h - hir ~ s) 

3:2 — xi , 16 
, In 

y'XlX2 



1 - k2 



We can also confirm the behaviour as r 
For r > s, we have 



(36) 



and r — + 00. 



£2 _ 

Xi 

and 



i-fi(i±i^4i^ + o(.-) 



3:2 xi _ 1 (4 + s")(4- 
Xi y X2 2 



0{r- 



(37) 



(38) 
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It follows from Eqn. @ that both (1 + Ai) and (1 + A2) go 
to 1 as r 00, so that the difference of Lambda functions 
in Eqn. ( ^5| ) goes to zero. The factor k ^ in this limit, 
and K( ^) ie finito, but tho oooflioiont of K gooe to aoro, from 
Eqn. (r 



sn /"i — In 



n anri A, ^ Oirr ga ovpoftorl 



For r < s, both xi and X2 diverge as r — ^ 0, but xijxx = 
1 + + O(r^), so that k ^ 0. Both K and Ao are finite 
here, so that I\ — I2 does not diverge, and the singularity is 
confined to the coordinates Xi. 

We may now move on to the difference Qi — Q2- After 
rewritin g If 2 in terms of K and Ao, the coefficient of K{hi) 
in Eqn. @ is {qk+q^v /(l + \2)) /2{xiX2f'^ ■ This is a rather 
messy expression in general, but at r = s, where x\ = 1, its 
value is 2{y/x2 — 1/ y/x2), so that the leading order term in 
Eqn. (H) is 



Aqjr s;s,c 
—r cos 2(j) 



2ir. 



1x2 - 



In- 



16 



1-k2 



(39) 



One can draw the same conclusion directly fro m Eqn. (|33[) 



by usin g the us eful asymptotic expansions in (Carlson and 



Gustaffon 1994), specifically relations (26), (34) and (44) in 
that paper. 

Since the singularity in Ai and Aq is no worse than 
logarithmic, we may numerically evaluate integrals involving 
these by using Gaussian quadrature with a log weight. 
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5 CONCLUSION 

I have obtained analytic angular integrals of the microlens- 
ing amplification function, for the case of a rotationally 
symmetric source. This avoids the need to use approximate 
methods to obtain this expression, and means that they can 
be evaluated more efficiently than using general numerical 
integrations. Also, we are able to make analytic statements 
about the leading-order behaviour of the integrals along 
their r = s singularity, and so use such asymptotic approxi- 
mations in further treatments. 

This also means that the dependence of the observed 
flux on the limb-darkening function, and of the observed 
polarization on the limb-polarization function, can be ex- 
pressed as integral equations. Thus the problem of recov- 
ering the latter from the former can be viewed as a classic 
inverse problem, which can be analysed in detail using the 
sophisticated techniques developed for such problems. This 
is the subject of a forthcoming paper (Gray and Coleman 



200c ) 



REFERENCES 

Abramowitz, M. and Stegun, I. A. (eds.), 1965, Handbook of 

Mathematical Functions, Dover, New York, Ninth printing 
Carlson, B. C, 1988, Mathematics of Computation 51(183), 267 
Carlson, B. C. and Gustafson, J. L., 1994, SIAM J. Math. Anal. 
25(2), 288 

Gaudi, B. S. and Gould, A., 1999, ApJ 513, 619, [astro-ph/9802205 
Gradshteyn, I. S. and Ryzhik, I. M., 1994, Tabic of Integrals, 

Series, and Products, Academic Press, 5th edition 
Gray, N., 2000, Automatic Reduction of Elliptic Integrals using 

Carlson's Relations, Submitted 



© 0000 RAS, MNRAS 000, 000-000 



